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Using the operator product expansion (OPE) technique, we study the rare double radiative decay 
B Kyy in the Standard Model (SM) and beyond. We estimate the short-distance (SD) contri- 
bution to the decay amplitude in a region of the phase space which is around the point where all 
decay products have energy ~ m/,/3 in the rest frame of the B-meson. At lowest order in l/m^, 
the B Kyy matrix element is then expressed in terms of the usual B ^ K form factors known 
from semileptonic rare decays. The integrated SD branching ratio in the SM in the OPE region 
turns out to be A.S§{B Kyy)^^^ ~ 1 x 10^^. We work out the di-photon invariant mass distri- 
bution with and without the resonant background through B K{ rjc , XcO } — * Kyy. In the SM, the 
resonance contribution is dominant in the region of phase space where the OPE is valid. On the 
other hand, the present experimental upper limit on B^ T+T^ decays, leaves considerable room 
for New Physics (NP) in the one-particle-irreducible contribution to B ^ Kyy decays. In this 
case, we find that the SD B Kyy branching ratio can be enhanced by one order of magnitude 
with respect to its SM value and the SD contribution can lie outside of the resonance peaks. 
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Flavor changing neutral cuiTents (FCNC) S-decays are of utmost importance in exploring and 
shaping the flavor structure of the SM and its extensions, since they ai^e forbidden in the Bom 
approximation and highly sensitive to NP contributions in the loops. Although many studies, ei- 
ther experimentally or theoretically, have been devoted to — > 57 and b —>■ si^i^ transitions, less 
attention has been paid to the double-photon decays mediated by — > syj. At the quark level it 
receives one-particle-reducible (IPR) contributions from the b ^ sy transition plus an additional 
photon and a one-particle-irreducible (IPI) term from a fermion loop with the two photons emitted 
from that loop (see Fig|T]). At the lowest order, the corresponding amplitude can be written as [1,2] 

■Uh{p)e^{ki)Ev{k2), (1) 

where A; = VthV^* and p{p') represents the momentum of the ft(5')-quark. We denote by k\,k2 
the 4-momenta and by £(^1,2) the polarization vectors of the photons. The tensors w!^^ and W^^ 
exhibit the contribution of the four-fermion operators ^1,2 — {sY^Lc){cY^Lb) and the photon dipole 
operator 1^7 ~ sG^yF^^Rb, respectively, and are given together with the loop function Kc in [1]. 

In this work we consider the exclusive B —>■ Kyy decay. Since its matrix element induced by 
^7 is non-local, it is difficult to estimate model-independently. Previous works [3,4] used here 
solely vector meson dominance models. Instead, we use OPE techniques, which allow to / make 
the matrix element local in a specific kinematic region by choosing appropriate photon energies 
and then // express the matrix element in terms of known form factors. Specifically, we demand 
both the internal s- and Z^-quarks to be far off-shell, with virtualities of order nih- For example, if 
both photons have energies nih/S, then the intermediate propagators of the IPR diagrams {Q\ 2)^ = 
{p' + )^ and (2j 2)^ — ~iiP ~ ^i^^)^ ~ "^fo) equal mi,/ y/3 and mh^Jlfi, respectively. Then we 
integrate out large scales of order m^. We construct the vertices in the effective theory out of a 
bottom heavy HQET quark /i,, and a strange coUinear SCET quark % [5,6]. Here, v = psh^B and a 
light-like vector n = pk/Ek, where Ek,Pk denotes the energy, 4-momentum of the kaon. Hence, 
we perform an OPE in Aqcd / Q, where Q = {mh,EK, Qt\2^ \f?} and = {k\ + k2)^- 

For the lowest order matching onto the b syy amplitude, we obtain after using the equations 
of motion the following dimension 8 operators [1] 

'^XO^vC>aiiRh,FfF^\ -2imao^,yRh,FrFfv''nP, U^Lh,FfDaF2fi^ +(1^2) 

where F^^ = 1 /2£^va/3^"^ and Wilson lines are understood in %■ The B Kyy matrix elements 
are then obtained from tree level matching of the QCD onto the SCET currents [5,7,8] as 

{K{n)\xh„\B{v))=2EK^iEK), {K{n)\xo^yh,\B{v)) = -2iEKi;{EK) (v^«v -Vv?v) , 

{K{n)\xy^h\B{v)) = lEKi;{EK)n^, {K{n)\xo^yYsK,\B{v)) = -2EKaEK)e^vai}v''nP , (2) 

where the form factor ^{Ek) can be identified with the QCD form factor in the usual parametriza- 
tion, as = /+. It is the only form factor remaining in the symmetry limit, known e.g. from studies 
in B ^ Ki^£^ decays. Analytical formulae for the B Kyy amplitude from the OPE are given 
in [1]. 

An additional contribution to B ^ Kyy beyond the OPE is photon radiation off the spectator. 
It is either kinematicaUy excluded (soft gluons between active quarks and spectator) or suppressed 
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Figure 1: Leading order Feynman diagrams for b — *■ syj decays. Left plot: IPI diagram from four-fermion 
operators, other two plots: insertion of O-]. Figure taken from [1]. 



by ttj. Annihilation contributions are leading power in l/Q, but suppressed by CKM elements or 
small penguin Wilson coefficients (for further details, see [1]). 

Long-distance effects from B (t]v YY)K, where T]v = T], T]', T]c and B {K* Ky)y are 
sizeable in B ^ Kyy decays [3,4]. Both the K*- and the T]^') -channels require kinematics that is 
outside of our OPE window. Therefore, we take in our analysis only the charmonia, specifically 
the r\c and also the XcQ^ into account. (The current upper bound on .'^{B —>■ XciK) implies that the 
background from the Xc2 is at least one order of magnitude down with respect to the scalar one.) We 
estimate the long-distance contamination from rjc, XcO with a Breit-Wigner ansatz fitted to data [1]. 

The resulting B Kyy di-photon mass spectra dT/dq^ are shown in Fig. ^ (For cuts and 
input, see [1]). The SD part of the SM spectrum (solid curve) is completely hidden behind the 
Tjc and the XcO resonance contribution (dash-dotted curve). The integrated SD branching ratio in 
the SM with cuts is small, A^(B Kyy)^^^ ~ 1 x 10"^, with about 50 % uncertainty from the 
renormalization scale. The contribution from IPI SD diagrams only gives a somewhat reduced 
branching ratio Kyy)^^^^^' ~ 0.5 x 10^^. Such small SD branching ratios for double 

radiative decays are not a complete surprise. In fact, being of the same order in a^m as semileptonic 
rare decays with branching ratios of few x 10^^ [9], we expect with C2 ^ 1.1 and Cf^ ~ —0.3 
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The double photon decays are substantially suppressed with respect to the semileptonic ones, since 
the di-lepton operators with large coefficients |Cg^| ~ ^(4 — 5) ai^e not contributing to the former. 

To be compai^able to the resonance contributions at least in some region of phase space the SD 
contribution has to be lifted by roughly one order of magnitude above the SM one. At the same 
time, such type of NP should not violate other data. A sizeable enhancement of the photon dipole 
coupling, C7, is excluded by data on B ^ X^y, which forces \C^\ ~ \Cj^\. More room for NP is in 
the 4-Fermi operators, which contribute at higher order to ft —> 57 decays. A possibility are large 
non-standard (pseudo)scalar couplings to taus ^^^^p^ = ^sRM{y=,)x. From =^(5^ T+T") < 5% 
[10] one obtains jCj^^^ | < 700, which is also consistent with data on other rare decays such as ft —> 57 
and ft — > sl'^l^ , 1 = e,}JL. In particular, large couplings ai^e allowed because FCNC decays into a 
tau pair are essentially unconstrained to date. In other words, decays such as B^ —>■ T+T^ and B — > 
K^*^T^T^ have sizable room for NP. As can be seen from Fig.^, the corresponding maximal SD 
B Kyy spectrum with NP in the T-couplings leaks mai^ginally outside the resonance background. 

Having a close look at the exclusive B —>■ Kyy decays, using the HQET and the SCET for- 
malisms, we found that the resulting SD branching ratio in the SM is small, order 10^^, in agree- 
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Figure 2: Di-photon mass distribution of B ^ Kyy decays in the OPE region. The solid (dashed) curve 
denote the SM pure SD contribution with ( and without) the IPR terms. The dotted ( double-dotted) curve 
corresponds to a NP scenario with CJ, = — 700(+700) andC^ = 0. The dash-dotted curve is the SM including 
the resonance contributions from the rjc and the XcQ- Figure taken from [1 ]. 



ment with the expectations based on the semileptonic rare B K£^£^ decays. We discussed 
further contributions, including long-distance effects via B — > rj^K — > jjK and B — > XcoK JYK. 
Unfortunately, we found that the SM SD contribution is not accessible behind the resonance peaks. 
We have further explored the possibility of NP contributions in B — > Kjj. In particular, enhanced 
FCNC couplings to taus can give via the IPI-loop branching ratios up to ~ 10^^. The correspond- 
ing maximal SD di-photon spectrum is then for some range of above the resonance background. 
We conclude that it is difficult to test SD physics with exclusive B Kyy decays. However, even 
strongly resonance polluted, this mode could be useful to probe the photon helicity in — > 57 via 
interference of different resonance amplitudes, as pointed out recently in [11]. 
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